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ABSTRACT 

We point out that the usual self-similarity in cold dark matter models is 
broken by encounters with individual normal galactic stars on sub-pc scale. Tidal 
heating and stripping must have redefined the density and velocity structures of 
the population of the Earth-mass dark matter halos, which are likely to have been 
the first bound structures to form in the Universe. The disruption rate depends 
strongly on galaxy types and the orbital distribution of the microhalos; in the 
Milky Way, stochastic radial orbits are destroyed first by stars in the triaxial 
bulge, microhalos on non-planar retrograde orbits with large pericenters and/or 
apocenters survive the longest. The final microhalo distribution in the solar 
neighborhood is better described as a superposition of filamentry microstreams 
rather than as a set of discrete spherical clumps in an otherwise homogeneous 
medium. This has important consequences to our detections of microhalos by 
direct recoil signal and indirect annihilation signal. 

Subject headings: dark matter 



1. Introduction 

Cold dark matter (CDM) has had considerable success in accounting for the observed 
large-scale structure of our Universe. On galactic and sub-galactic scales, however, the abun- 
dance of low-mass structure and the degree of concentration of dark matter halos predicted by 
CDM simulations have provoked a considerable amount of discussion. To clarify how CDM 
behaves on the smallest scales, Diemand and collaborators recently performed the highest 
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resolution numerical simulations of dark matter clustering to date (Diemand et al. 2005). 
These simulations start with an initial spectrum of density fluctuations extending down to 
the free-streaming mass of the dark matter candidate, which for a generic weakly interacting 
100 GeV particle (such as a supersymmetric neutralino) is roughly ~ 10^® (t/o/10~®s)^/"^ Mq, 
where freeze-out occurs at an epoch t/o ~ 10^^ (m^/lOOGcV) s. As expected from scaling 
the results of simulations on larger scales, fluctuations in this initial distribution collapse and 
virialize to form structures with roughly 200 times the background density at a redshift of 
z ~ 50. This first generation of 'microhalos' survives to some degree as substructure in the 
larger halos that form subsequently in the simulations. Extrapolating to the present-day, 
Diemand et al. 2005 suggest that 10^^ Earth-mass clumps should survive in the halo of the 
Milky- Way, amounting to about 0.1% of its total mass. If this were the case, the nearest 
Earth- mass clump would be on the order of 0.1 pc distant from the Earth and would have a 
typical size of 0.01 pc. Such nearby clumps might be observable as sources of gamma-rays, 
produced in dark matter annihilations. They might also contribute to diffuse cosmic ray 
fluxes of positrons or anti-protons. In order to motivate experimental searches, however, it 
is important to determine whether the predicted small clumps of dark matter would actually 
survive the grainy tidal field due to stars while crossing the Galactic bulge and disk ~ 100 
times. 

Here we estimate the rate at which Earth-mass clumps are disrupted in stellar encoun- 
ters. We show that the nearest clumps in a galaxy hke the Milky Way are hkely to be 
tidally disrupted by repeated encounters over a Hubble time. Nevertheless, tidal debris from 
individual clumps may produce distinct microstreams that are potentially observable in dark 
matter detection experiments. 

2. A Semi-analytic Mass-loss Estimate 

Extended bodies generally have a complex response to external tidal heating (Gnedin & 
Ostriker 1999). We can estimate the net effect of this process by applying the basic scahng 
from the semi-analytic mass- loss model of Taylor & Babul (2001), which has been shown to 
match tidal heating and mass loss rates in high-resolution numerical simulations. According 
to this model, whenever a system on a general orbit spends time At in a strong tidal field, 
that is a field in which the instantaneous tidal limit for the object r^{t) (calculated as in 
the circular-orbit case - cf. Binney & Tremaine 1987) is smaller than its size r, it loses 
a fraction At/torb of the mass outside rt(t), where torb = 27ir{t)/v{t) is the instantaneous 
orbital period. 

To estimate when material at radius r in a microhalo will experience strong shocks, we 
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need to calculate the minimum impact parameter such that rt{t) < r. This is equivalent to 
the condition that the tidal force generated during an encounter with a star of mass at 
an impact parameter bsh, 2G'm*x/6^^, exceeds the restoration force for a small displacement 
X, ^^^^^x, where p^r is the mean density of the microhalo interior to r. For a typical 
microhalo in the simulations by Diemand et al. (2005), the density at the half-mass radius 
is ~ 1Mq])c~^, while the highest resolved density (at 0.1 r2oo) is ~ 10 times higher. We 
do not expect a divergent (phase space) density in these smallest halos. Even if a cusp 
existed below the simulation resolution, it would contain less than 3% of the microhalo's 
mass, extrapolating from the profile given by Diemand et al. (2005). Thus for a rapid 
encounter with a solar- mass star, impact parameters of bgh ~ {'itT'*/ Pmic)^^^ ~ 0.5 pc or less 
will produce a strong shock over a range of radii containing 97% or more of the mass in the 
microhalo. pmic is the average density within 10% of the virial radius -Rvir and can be thought 
of as a constant until the microhalo is virtually destroyed as shown in Angus & Zhao (2006, 
hereafter AZ, figures 7c and 12b). 

A microhalo moving at a relative speed of Vr ~ 300kms~^ through a star field with 
a number density = p^/m^ will encounter stars with an impact parameter < bgh at a 
rate of N^nc = n*l^7r6^^ ~ n=f/0.004pc^^ per million years. Each encounter shocks the 
microhalo for a time At = 2b/Vr ~ 3300 years, and liberates a fraction of the microhalo's 
total mass. This fraction, from the scahng in Taylor and Babul (2001), is roughly Am/m — 
At/toTh — forb/TrK ~ 10~^. Here, torb = 27765/1/ Vorb ~ 30 Myr is the period of a circular 
orbit of speed Vorh — -s/Gm^/bsh ~ 90m/s around the solar-mass perturber at a distance 
of bsh ~ 0.5 pc. Thus, the mass of a microhalo decreases exponentially, with the e-folding 
time given by 1/r = —m~^dm{t)/dt = NcncAm/m = Gn^m^{Gpiaic)~^^'^, where we have 
expressed b in terms of the parameters of the satellite. Averaged over a Hubble time, we 
find that a microhalo's mass decays exponentially, m{t) = moexp(— t/r), where the e-folding 
time and the present mass are given by 

1^ ^0 ^ lOGyr ^ f jp*) \ ( Pmic y'^' 

mnow T ~ VO-OO4M0pc-3y VlOMopc-V ' ^' 

where (p*) is the mass density of galactic (disk, bulge and halo) stars averaged along the 
orbit of a microhalo over the past 10 Gyrs (at earlier times, few stars might have formed). In 
general, (p*) is a function of the pericenter, apocenter and vertical height of the microhalo 
orbit. E.g., a microhalo on a thick disk orbit of height 1000 pc near the solar neighborhood 
would have made ~ 100 disk crossings and gone through star fields of an accumulated column 
density ~ 10'^ pc"^ by today. It goes through the galactic disk with a mean surface density 
of 46M0pc^^ or an average volume density 46M0pc^^/2OOO pc ~ 0.023MqPC~^, thus this 
microhalo's mass is predicted to have decayed by about 6 e-foldings. 
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Fig. 1. — shows the post-impulse fraction of bound microhalo particles vs. the equivalent 
column density S^, = M^/inh"^) of the perturbers, where symbols denote results from a set of 
numerical experiments of stellar encounters with microhalos of initially tangential anisotropy 
(TA, red symbols) nd radial anisotropy (RA, black symbols). In these experiments a per- 
turber with = O.IMq (plus), O.SMq (asterix), O.GMq (diamoNds) and IMq (triangles) 
passes a 10~^Mq microhalo with relative velocity = 200kms~^ and different values of 
the impact parameter b ranging from penetrative {b < Rm,- to non-penetrative encounters 
(6 > Rvir)- For the same mass column density, many low-mass penetrative perturbers are 
not as efficient as few high-mass non-penetrative perturbers in terms of stripping the mi- 
crohalo; stripping is unimportant for sub-stellar perturbers in general. Our semi-analytical 
prediction assuming pmic = lOM© pc~^ is represented as the solid line. 
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Another way to think of the semi-analytical disruption rate is to note that the logarithm 
of the bound fraction In "^"'^ oc f p^:dt oc y^, where is the stellar mass density encoun- 
tered. To get an indication of how robust the semi-analytical massloss rate is, we have also 
run N-body simulations of a simplistic single encounter between a star and a microhalo for 
comparison. We use the tree-code of Vine and Sigurdsson 1998. We use a total particle num- 
ber of 10^ and softening of 10~^pc consistent with interparticle separations of 3.5xl0~'^pc. 
The tree code tolerance parameter, 9 is the cell size divided by particle separation and is set 
to 0.6 in these simulations. Initially the microhalo particles were distributed in a spherical 
NFW density profile with a i?vir = 1.6Rc — O.Olpc and a total mass of 10~^Mq. Two sets 
of velocity space initial conditions were used; tangentially anisotropic (TA) microhalos with 
randomly inclined circular orbits V^{r) = and radially anisotropic (RA) microhalos with 

K^(^) — ^~2^ — 7^ X^""^ K^pc?^) where p{r) and K are the microhalo's internal density and 
circular velocity respectively (cf. Angus and Zhao 2006). Two-body relaxation time, mea- 
sured by the mean squared normalised energy drift after 10^ time steps using 10^ particles 
for an isolated microhalo, is about 3Gyr. So a microhalo resolved with 10^ particles should 
be free from numerical relaxation on time scales of 15Myr, i.e., the dynamical time scale of 
the microhalo. To measure the effect of encounters with stars, we shoot a single star of mass 
M* = 0.1, 0.3, 0.6, IM© with a relative velocity Vr — 200kms~^ towards a microhalo with 
a closest approach & = (0.2 — 1) x i?vir for penetrative encounters And 6 = (1 — 4) x i?vir 
for non-penetrative encounters. Such an impact is typical for a microhalo passing through a 
random star field with an equivalent star column density = M*/(7r6^). The duration of a 
star passage, say, along the x-axis from x = 10i?vir to x = — 10-Rvir is about 0.0001 dynamical 
time of the microhalo, hence microhalo particles are virtually stationary (moving less than 
0.01 percent of its orbital radius) during a star passage, so the perturbations are indeed in 
the impulse regime. It appears that our semi- analytical model captures the essential feature 
in numerical simulations of the impulse on 10~^Mq microhalos if assuming pmic ~ l-^o pc~^ 
(cf. Fig.l). However, there are microhalos of a few times Earth mass, and these can be, e.g., 

5 times more massive and denser than the 10~^ Mq microhalos experimented here (cf. Fig. 2 
of Diemand et al. 2005). Allowing for more robust microhalos and to be conservative in 
estimating the disruption rate of microhalos, we shall adopt the rate from the semi-analytical 
model assuming pmic ~ IOMq pc~^ hereafter. Note that the scaling of our semi-analytical 
model Eq. 1 also agrees with the literature on impulse-induced massloss in the context of 
diffuse star clusters, e.g., Moore (1996), and Eqs. 6-8 of Wielen (1985) or Eq. 8 of Wasserman 

6 Salpeter (1994). 
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Fig. 2. — Left panel: The inset shows the scatter diagram of the velocity space (at launch) 
for the 500 microhalos launched from the solar circle. Microhalos are given different symbols 
depending on the time-averaged stellar density along the orbits ((p*) > O.OIMqPc"^ for 
crosses and < O.OOIMqPc"'^ for solid circles and in between for squares). Note microhalos 
with low orbital energy and angular momentum are first to be disrupted. Three typical 
orbits (larger symbols) are picked out and shown in meridional plane Z vs. R. The yellow 
shaded area is a cut of the Besancon star count model of the Galaxy (Robin et al. 2003) 
with equal star density of 10~*Mq pc~^ for i = 0,1,2,3 plotted in contours. Right Panel- 
shows the star density in Mq pc~^ along the three colored orbits as function of time. 
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3. Microhalos on Near-Ecirth Orbits 

To study the dependence of mass loss on orbital parameters, we have integrated orbits 
of microhalos in a flattened, axisymmetric galaxy potential with a nearly flat rotation curve: 

^{R,z) = (220kms~ )^ In ^^:^^..i^gQ^,pg - We launch 500 orbits from the solar neighborhood, 
{R, z) = (8kpc,0), with a Gaussian velocity distribution (150 kms^^ dispersion in each 
direction). We read off the stellar density along the orbit using the Besancon model of the 
Galaxy (Robin et al. 2003), and then take a time average of the star density. The Besancon 
model is a widely accepted detailed prescription of the number density of stars and stellar 
remnants of different ages and masses in the Galactic spheroids (bulge and halo) and disks 
(thin and thick). A few examples of the microhalo orbits are given together with the star 
count model in Fig. 2. These simulations suggest that the population of microhalos with 
apocenters less than 12 kpc has (p*) ~ O.OO2M0pc^^, and are hence past their half-life 
(purple orbit in Fig. 2). The more planar orbits, with \Z\ < 4(/?/20)kpc and apocenter 
R < 20 kpc (blue orbit in Fig. 2), have decayed by approximately 1 e-fold with (p^,) ~ 
O.OO4M0pc"^. Complete destruction should happen to disky orbits penetrating into the 
bulge; e.g., the red orbit in Fig. 2 goes through a typical column density of SOOM© pc"^ with 
each disk crossing at a small angle, and has on average (p*)^ 0.075Mq pc~^, corresponding 
to about 19 e-folds. In general, we observe a strong correlation between the orbital shape 
and the disruption rate. Orbits which go through dense regions of the Galaxy are likely 
disrupted due to tides (cf. Fig. 2). 

Also plotted in the inset of Fig. 2 is the zone of destruction in the velocity space using 500 
microhalos launched from the solar neighborhood. There is a clear trend for halos with more 
planar (smaller \Vz\) orbits, especially those co- rotating with the disk, and orbits with lower 
energy (hence smaller apocenters) to enter regions of higher density. The sample averaged 
density is O.OOSMq pc~^ or 0.02Mq pc~^ for median or mean density, so microhalos on Earth- 
crossing orbits are likely to have been severely stripped. About 10% of the microhalos are 
in low-density regions (below O.OOIM© pc~^), and might survive more-or-less intact. The 
chance of survival, however, is likely even more severely reduced in the triaxial, barred and 
evolving potential of the Milky Way, where most orbits are stochastic box orbits which pass 
through the dense center of the Milky Way at some point over a Hubble time. Additionally, 
it has been shown in AZ that disk shocking and tidal stripping may combine with stellar 
encounters to destroy microhalos even more efficiently than shown here. 

Of course, these orbits are mearly a sample of all orbits passing the sun. In fact, the 
fraction of halos disrupted would be much lower for high-energy orbits beyond the solar 
neighbourhood as the tides at larger pericentre are weaker along with the rapidly decreasing 
fraction of stars. What was also keenly noted by Berezinsky et al. (2005) was that apocentric 
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disk crossings beyond the extent of the thin disk are hkely to boost the hfetime of microhalos. 
The destruction rate is certainly much higher for low-energy orbits contained inside the bulge 
or those that have pericentres within the solar radius. These michohalos would be exposed 
to stronger tides and higher stellar densities added to the fact they will cross the disk 
twice every orbit. However, these orbits are irrelavent for direct experiments on Earth, and 
less important for indirect experiments because microhalos are barely identifiable from the 
background beyond about O.lpc - Ipc distance from us. 

Additionally, the calculations done by Diemand et al. (2005) did not find self-consistent 
orbits and their comments on the current situation in the Milky Way are major extrapo- 
lations from dark- matter only simulations that stop at redshift 26, and only ever simulate 
a much smaller total volume. Since halo orbits in the MW would be strongly affected by 
the formation of the galaxy. So what we have here is desruption as a function of orbital 
parameters. 

We note that the long term fate of the cusp at the very center of each microhalo is 
unclear. Simulations have shown that systems with a universal density profile may lose all 
but 0.1-1% of their mass and still retain a bound central region Hayashi et al. (2003). 
Such objects may indeed survive to the present-day in the solar neighborhood, but they will 
have little effect on dark matter detection experiments. The tidal debris stripped out of 
microhalos may actually be of greater interest, as discussed below. Nevertheless, AZ (figures 
7c and 12b) showed that the inner density of a microhalo (pmic) is largely unaffected until it 
is stripped down to the inner shell. 

4. Implications For Direct and Indirect Detection 

The existence of dark substructure in our galaxy's halo has important implications for 
the prospects of both direct and indirect dark matter detection (Green et al. 2005; Berczinsky 
et al. 2003; Fieri et al. 2005; Koushiappas et al. 2004). The presence of substructure 
affects direct detection rates in a straight forward way: If the solar system happens to be 
located inside of an overdense region of dark matter, then the rate observed will be enhanced 
proportionally to the density of the clump. If, as is much more likely, our solar system is not 
inside a clump, then the rate will be modestly reduced by the fact that some fraction of the 
overall density is contained in substructures and thus does not contribute to the density of 
the smooth component. 

Zhao and Silk 2005 estimated that the nearest minihalo of lO^M© is 1 kpc away with a 
size of 50pc. Diemand et al. 2005 estimated that unstripped Earth-mass halos have a filling 
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factor of 0.5%. Hence in either case it is unlikely for the Sun to be inside a substructure. 

Tidal debris from encounters with stars changes the picture substantially. While the 
probability to find our solar system inside of a lO^^pc-sized bound remnant of a microhalo 
is on the order of 10~^, the probability is of order unity for the solar system to be moving 
presently inside an (unbound) substructure. Disrupted streams are inhomogeneous regions 
with at least 10 times the volume filling factor of the original unstripped clump. This is 
because at each disk crossing a burst of tidal debris is released due to encounters with disk 
stars, and the tails would lead/trail at least O.lpc (assuming Im/s speed of escape) along 
the system's trajectory in 100 Myrs, the time between disk crossings. The time for the 
solar system to cross a stream at a general angle will be longer than 50 years, the time to 
cross a bound microhalo (Diemand et al. 2005). Moving inside the inhomogeneous debris 
of the microhalos could cause transient enhancements of the event rate in direct detection 
experiments. Furthermore clumps on polar orbits survive longer than planar orbits (cf. 
Fig. 2), this makes the velocity distribution of the micro-streams non-Gaussian, and could 
imprint an interesting annual modulation on the direct detection signal (Morgan et al. 2005; 
Freese et al. 2005). Fig. 12 of Angus & Zhao (2006) shows the morphology of a microstream: 
even after the microhalo is fully disrupted, the microstream is still a distinct entity with a 
negligable tangential size compared to its length. As a result, these filamentary streams may 
be the fate of a large proportion of the Milky Way's microhalos and the enhanced background 
from these structures may impose interesting features on direct detection experiments as well 
as having interesting annihilation flux properties. 

The implications of dark microhalos and microstreams are rather different for the case 
of indirect detection. Indirect measurements sample the distribution of dark matter, through 
its annihilation rate, over large regions of the halo. Furthermore, dark matter annihilation 
rates, and thus indirect detection signals, are proportional not to the density of dark matter 
but to the dark matter density squared. Substructure within the galactic halo thus may 
be capable of boosting the dark matter annihilation rate and enhancing the prospects of 
detecting dark matter indirectly (Taylor & Silk 2003). 

Techniques employed for the indirect detection of dark matter include gamma-ray, anti- 
matter and neutrino detectors (Silk & Srednicki 1984). Gamma-rays annihilating in nearby 
and dense substructures, if present, could provide point sources potentially observable by 
gamma-ray telescopes (Fieri et al. 2005; Koushiappas et al. 2004). Anti-matter produced in 
dark matter annihilations are deflected by galactic magnetic flelds and thus only the diffuse 
spectrum can be studied. Neutrinos are not as useful for identifying annihilations in the 
galactic halo, but instead are used to search for dark matter particles annihilating in the 
core of the Sun. The neutrino flux produced through dark matter matter annihilations in 



-10- 



the Sun is tied to the local density averaged over very long periods of time and therefore the 
presence of substructure is of little importance. 

Nearby microhalos, being tidally disrupted, are not particularly bright point sources 
of gamma-rays. While Diemand et al. 2005 predict that the nearest intact microhalo will 
be much brighter than known dwarf spheroidals in gamma-rays, we find that the situation 
reverses for a microhalo left with only a one percent bound remnant. 

Anti-matter fluxes produced through dark matter annihilations do not depend on the 
nature or location of individual dark substructures, but rather on the distribution of dark 
matter averaged over large volumes. The value of the quantity (p^) / (p)^, averaged over the 
contributing volume (a few kiloparces for positrons or tens of kiloparsecs for anti-protons), 
determines the overall boost factor for the anti-matter fluxes generated through dark matter 
annihilations. 

With a distribution of substructure of the form dN/d\ogM oc M~^, each decade of 
mass contributes almost equally to the boost factor, so it is important to determine over 
what range of masses substructures can survive. Previous works (Helmi et al. 2002; Stochr 
et al. 2003) emphasized the role of largers substructures, with the mass of Draco dwarf 
spheroidal. Diemand et al. 2005 argue that the nearest Earth-mass halo is factor of a few 
brighter than Draco. We argue, however, that the nearest Earth-mass halo is likely to be 
tidally dispersed, and hence fainter than Draco. We have proposed (Zhao & Silk 2005) that 
the nurseries of the first stars in the universe are dense minihalos of 10^ M© formed at 2; ~ 22. 
Unlike the microhalos such minihalos can survive the galactic tides until the present day. 
The nearest minihalo is only 1-2 kpc from the Sun, and is an order of magnitude brighter 
than the combined annihilation from the Draco galaxy. 

In hght of the excess reported by the HEAT collaboration (Barwick et al. 1997; Coutu 
et al. 1999; Baltz & Edsjo 1999), substructure is of particular importance for dark matter 
searches involving positrons. To produce such a signal with thermally generated neutralinos 
or other type of WIMPs, however, boost factors of ~50 or higher are required (Baltz et al. 
2005; Hooper & Kribs 2004; Hooper & Servant 2005). Although it has been shown that 
more massive substructures do not naturally provide such a large boost factor (Hooper et 
al. 2004), effects due to baryons, or the presence of intermediate mass black holes could 
perhaps make them brighter (Prada et al. 2004), although probably not enough to account 
for the HEAT excess. After the effects of tidal disruption are considered, the same is found 
to be true for Earth-mass microhalos. We agree with previous studies (Berezinsky et al. 
2003) which determine that positron boost factors larger than ~2-5 are unlikely. Given this 
conclusion, another explanation for the excess observed by HEAT appears to be required. 
Our understanding of these issues will be dramatically improved with data from the up- 
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coming cosmic anti-matter experiments PAMELA and AMS-02 (Profumo & UUio 2004). 

5. Discussion and Conclusions 

Using a semi-analytic estimate of the tidal mass-loss rate, together with numerical in- 
tegration of a fair sample of microhalo orbits, we find that most microhalos present in the 
solar neighborhood will have been heavily stripped by stellar encounters, producing 'mi- 
crostreams' of tidal debris. More generally, in environments with very low stellar densities 
such as the outer parts of the Galactic disk or in a Scxtans-like dwarf galaxy, microhalos 
are only mildly heated. But in high-density environments such as the Galactic bulge or an 
M32-like elliptical (800 to O.OSMopc"^ in inner Ikpc of M32 (Mateo 1998)) microhalos are 
likely fully destroyed. In disk galaxies such as the Milky Way, the fraction of microhalos 
which are disrupted depends strongly on the orbital inclination and pericentre. Microhalos 
on orbits coplanar with the disk are very quickly disrupted. Microhalos in the outer part 
of galaxy halos will encounter far fewer stars, however, and thus their annihilation rates 
will be less dramatically reduced. A possible consequence of this might be to detect little 
gamma-ray emission from the central part of external galaxies, where stars have disrupted 
most of the dark substructure, and more flux from the outer regions. A ring of gamma-rays 
surrounding a galaxy, if detected, would provide a strong confirmation of the existence of 
dark substructure. 

The tidal effects discussed in this paper and in several preprints appeared during the 
marathon referee process of this paper (Berezinsky, Dokuchacv, Eroshenko 2005, astro- 
ph/0511494, Green & Goodwin 2006, astro-ph/0604142, Gocrdt, Gncdin, Moore, Diemand, 
Stadcl 2006, astro-ph/0608495, Angus & Zhao 2006, astro-ph/0608580) have important con- 
sequences for the direct and indirect detection of dark matter. In particular, the probability 
of the Earth intersecting a tidal stream at any given time is considerably larger than the 
probability of it being inside a microhalo 's bound remnant, and thus tidal streams may po- 
tentially increase the reach of direct detection experiments. Indirect detection rates are not 
enhanced, however, as once stellar encounters are considered. Earth mass halos are not a 
significant contributor to the overall dark matter annihilation rate. 
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